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We report on the results of our calculation of meson correlators in a finite volume. The calculation 
is carried out in the quenched approximation near the chiral limit (down to m q = 2.6 MeV) using 
the overlap fermion. For these small quark masses, the scalar and pseudo-scalar correlators are 
well approximated with a few hundred eigenmodes. The results for both connected and discon- 
nected correlators are compared with the theoretical predictions of quenched chiral perturbation 
theory. 
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1. Introduction 



In a finite volume the meson correlators are largely distorted. Under an extreme condition such 
as the £ -regime [j|], the correlators do not even show the familiar exponential fall-off as a function 
of t. It is, however, still possible to analytically describe the system using the chiral perturbation 
theory (ChPT), as long as the pion mass is small enough. One-loop calculations in the quenched 
approximation are available for both connected and disconnected correlators [||, ^J. Since they are 
expressed in terms of the usual low energy constants (LECs) in ChPT, the lattice QCD may be used 
to determine their values through the calculation of meson correlators in the e-regime. 

In this work we calculate the quenched meson correlators near the chiral limit. We employ the 
overlap-Dirac operator and calculate low-lying eigenvalues and eigenvectors in order to investigate 
the low-lying mode saturation of the meson correlators. For the scalar and pseudo-scalar correlators 
the saturation is well satisfied, and we can obtain them very close to the chiral limit without extra 
computing cost. We present our result for the pion decay constant F n and chiral condensate E. A 
full detail of this work is given in [Q] (see also [§]). Other works in this direction include [§, [7[ ^J. 



2. Simulation methods 



We work on a 10 3 x 20 lattice at j8 = 5.85 in the quenched approximation. The physical lattice 
spacing corresponds to 0.123 fm. We use the overlap-Dirac operator 

D m ={\-^jD + m, D=i(l + y 5 sgn(a# w )), H w = y 5 \D w -^\, (2.1) 



with a = a/{\ + s). The parameter s is set to 0.6. In the computation of the sign function (2.1), we 
treat 60 lowest eigenmodes of Hw exactly and approximate the rest by the Chebyshev polynomial 
(degree 100-200). The topological sectors are identified unambiguously by counting the number of 
zero mode of D. We analyzed 20, 45, 44 and 24 gauge configurations for topological charge \Q\ = 
0, 1, 2 and 3, respectively. The 200 lowest eigenvalues and associated eigenvectors of the overlap- 
Dirac operator D are calculated using the implicitly restarted Arnoldi method (implemented in the 
ARPACK numerical package). 

Near the massless limit, the quark propagator is expected to be well approximated by these 
low-lying eigenvectors. 

200+ 1 ei ^ 

°?<*y)~ I ( i- a »/2)^ + „ , 'W , ?M- <2 - 2) 

where A,- and v/(x) are the i-th eigenvalue and eigenvector, respectively. In Figure p] we compare the 
pseudo-scalar (left panel) and axial-vector (right panel) correlators constructed from the approxi- 



mated quark propagator ( J2.2[ ) with those from the exact propagator (obtained using the conventional 
CG inverter). At a small quark mass am = 0.008, we find a good agreement of the approximate 
correlator with the exact one for the pseudo-scalar correlator. Such saturation becomes even better 
for smaller quark masses and higher topological sectors. The saturation is, on the other hand, much 
worse for the axial-vector correlator. Therefore, in the following analysis the axial-vector correla- 
tors are calculated with the exact method, while the low-mode approximation is used for quantities 
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Figure 1: Comparison of the pseudo-scalar (left) and axial- vector (right) correlators evaluated with the 
low-mode approximation (stars) with the exactly solved correlators (open squares). 

for which the saturation is well satisfied. One advantage of the low-mode approximation ( |2.2[ ) is 
that one can reduce the statistical noise drastically using the low-mode averaging technique [[7], [H]] 
without much extra computing cost. 

3. Axial- vector current correlator 

For the local axial-vector current A" = YYsY^i /2) W we calculate the renormalization factor 
Za non-perturbatively using the Ward-Takahashi identity. >From the ratio of (VoAq(Y).P 3 (0)) and 
(P 3 (t)P 3 (0)} correlators, Za is obtained as a linear slope as a function of m. Our result is Za = 
1.439(15). 

In Figure § we show the results for the axial- vector correlators for topological charges \Q\ = 
0, 1 and 2. The expected form from the quenched ChPT (QChPT) is written as 

Z 2 A (Al(t)A^) = I f+2mZ m (n)Th 1 (t/T). (3.1) 

where £g(jU) is the scalar condensate at a given topological charge Q and /I = mLV: £q(/i)/E = 
H(Iq(IJ.)Kq(ij.) + Qj\l. (/q(jU) and Kq{\i) denote the modified Bessel func- 

tions.) The t dependence is simply given by a quadratic function h\ (x) = 1/2 ( (x — 1 /2) 2 — 1/ 12) . 
>From this expression, F % is determined through the constant term, while the second term has de- 
pendence on \Q\ and \i as well as t/T. >From a simultaneous fit for the data at \Q\ = and 1 with 
several values of m and t/T, we obtain F % = 98.3(8.3) MeV and Z 1 / 3 = 259(50) MeV. The formula 
( [Q| ) explains the data very well at \Q\ = and 1, but the same set of parameters does not describe 
the data at higher topological charge (\Q\ =2, shown in the right panel of Figure This is not 
inconsistent with the theoretical expectation, because the derivation of ( |3.1[ ) assumes a condition 
\Q\ < (2 2 ) and for our small lattice (Q 2 ) = 4.34(22). For more detail, see Ml. 

4. Connected scalar and pseudo-scalar correlators 

The connected scalar and pseudo-scalar correlators are obtained using the low-mode approxi- 
mation. Since we can average the source point over the space-time lattice, the signal is very clean 
as shown in Figure |3[ 
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Figure 2: Axial-vector current correlator at Q = (top left), \Q\ = 1 (top right) and \Q\ = 2 (bottom). The 



lines show the fit with the function (3.1) for \Q\ = and 1 



The fitting with the QChPT expectation is more involved than the case of the axial-vector 
current correlator, because the parameters rn\ and a appear to describe the quenched artifacts. To 
avoid a four-parameter fit with (F n , E, wiq, a), we input F n from the result of the axial-vector current 
correlator and use the relation (Q 2 ) /V = m^F K /2N c to determine tHq. 

The data shown in Figure || are well described by the QChPT expectation for \Q\ = and 1, as 
in the axial-vector current correlator. We obtain I 1 / 3 = 257(14)(00) MeV and a = -4.5(1. 2)(0.2), 
where the first error is statistical while the second represents the uncertainty of input parameters. 

Our result for a is rather large and negative, which contradict with the small a observed in 
many previous lattice calculations (for example, see [11]). We also attempted a fit with F n and £ 
as free parameters, while a is set equal to zero. The result is F n - 136.9(5. 3)(0.9) MeV and Z 1 / 3 = 
250(13)(00) MeV; F % is not consistent with the analysis of the axial- vector current correlators. 



5. Disconnected pseudo-scalar correlator 

We also study the disconnected pseudo-scalar correlator using the low-mode approximation. 
For this quantity we did not check the consistency with the exact solution because it is too costly. 
Instead, we looked at the saturation as the number of eigenmodes included in the approximation. 
Figure [| shows that the disconnected diagram is well saturated with the number of eigenvalues 
greater than 10. 

Data of the disconnected pseudo-scalar correlators are shown in Figure |5[ As in the connected 
correlator, we fit to the QChPT formula with £ and a as free parameters. The result for \Q\ = and 



133/4 



Meson correlators in a finite volume near the chiral limit 



Kenji Ogawa 



lattice Q=0 
lattice Q=1 
lattice Q=2 
qChPT Q=0 
qChPT Q=1 
qChPT Q=2 



h m m m a 



10 
t 



12 14 16 16 



20 





lanice \J= i i i i 




lattice Q=2 >--«.»! 




qChPTQ=1 - 




qChPTQ=2 


k-. 





















10 12 14 16 18 20 
t 



Figure 3: Connected pseudo-scalar (left) and scalar (right ) correlators. The curves show the fitting with the 
analytically predicted form. 
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Figure 4: Saturation of the disconnected correlator at t = 7, 9 and 10 as a function of the number of 
eigenvalues in the lowmode approximation. 



1 data is I 1 / 3 = 227(32)(00) MeV and a = -3.5(1.2)(0.3), which is consistent with the connected 
correlators. The analysis with a = yields F % = 125.7(5.6)(0.9) MeV and I 1 / 3 = 223(29)(00) MeV. 



6. Summary 

We performed a quenched simulation in the £-regime at small quark masses m = 2.6-13 MeV. 
In this small quark mass region, the connected scalar and pseudo-sealer correlators are approxi- 
mated by the low-lying modes of the Dirac operator to 1 % accuracy or even better. In the small 
topological sector, the results are well described by the one-loop QChPT formula; by fitting we can 
extract F % and £ as well as ;mq and a. 

We also find problems: Data at higher topological sectors are not consistent with the QChPT 
expectation. Even for the small \Q\ sector, the large negative value of a may suggest a breakdown 
of the theoretical expression. We suspect that these problems arises from the small physical volume 
of our lattice for which the condition \Q\ <C {Q 2 } is not well satisfied except for <2 = 0- 
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Figure 5: Disconnected pseudo-scalar correlators for \Q\ = l(left) and \Q\ = 2(right). The lines are the 
fitting function with the parameters which was obtained by calculating axial-vector current correlators and 
connected pseudo-scalar and scalar correlators. 
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